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We investigate freezing characteristics of a finite dipolar hexagonal plate by the 
Monte Carlo simulation. The hexagonal plate is cut out from a piled triangular lattice 
of three layers with FCC-like (ABCABC) stacking structure. In the present study an 
annealing simulation is performed for the dipolar plate uniaxially compressed in the 
direction of layer-piling. We find spin melting and refreezing driven by the uniaxial 
compression. Each of the melting and refreezing corresponds one-to-one with a 
change of the ground states induced by compression. The freezing temperatures of 
the ground-state orders differ significantly from each other, which gives rise to the 
spin melting and refreezing of the present interest. We argue that these phenomena 
are originated by a finite size effect combined with peculiar anisotropic nature of the 
dipole-dipole interaction. 



I. INTRODUCTION 

In recent years, systems consisting of arrayed single-domain ferromagnetic nanoparticles 
receive much attention of many researchers because of its possibility as a high storage den- 
sity 0]. In such array systems the dipole-dipole interaction is considered to play an important 
role on their magnetic properties. The interaction not only intrinsically involves frustration, 
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i.e., the competition between ferromagnetic and antiferromagnetic interactions, but also it 
connects the anisotropy in the the spin space with that in the real space. In finite systems 
the latter is expected to yield a vast variety of ground-state spin configurations depending 
on shapes of the systems. Various ordered spin patterns have been actually observed, for 
example, in infinite P, 0, Ol or finite 0, 0] systems only with the dipole-dipole interaction, 
which we call dipolar systems, on the triangular lattice. In dipolar systems such a peculiar 
phenomenon as the from-Edge-to-Interior freezing has been also found recently In the 
present study we investigate a dipolar hexagonal plate cut out from a piled triangular lattice 
of three layers with FCC-like structures and with uniaxial compression in the direction of 
layer-piling. Our aim here is to demonstrate that, in the dipolar system even with such a 
restricted shape, an interesting crossover is predicted to exist when the interlayer distance / 
is changed. 

Figure HI shows a schematic spin-pattern and T* diagram obtained by the present study of 
the dipolar hexagonal plate, where T* is the temperature at which spins start to be frozen. 
At low temperatures there appear three kinds of patterns, i.e., in-plane- vortex, convex- flow, 
in-plane-af-FMC (see Fig. Ol^c) below) patterns depending on a value of /. The freezing 
temperature, T*, depends on a spin pattern of the corresponding ground state. Especially, 
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FIG. 1: The schematic spin-pattern and T* diagram of the dipolar hexagonal plate. Arrows in 
insets shows directions of spins. 
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T* of the convex flow state is much lower than those of other states. Because of this l- 
dependence of T*, we expect the peculiar spin melting and refreezing driven by uniaxial 
compression, say, at T = 0.3. 

In the present study we perform simulated annealing of the dipolar hexagonal plate by 
the Monte Carlo (MC) method. The well-known Hamiltonian of a dipolar system. 



is employed. Here S'j's denote classical Heisenberg spins which represent magnetic moments 
of nanoparticles, D a coupling constant, and a a lattice constant of the triangular lattice. 
The number of piled layers in a plate examined, Ni^, is 3. Each layer contains about 100 
sites. The ^-direction is set in perpendicular to the layers. We perform typically 8 cooling 
runs with a fixed cooling rate, in which we use different sets of random numbers generating 
their initial configurations as well as spin-renewing for each MC step. For each cooling run, 
temperature T is initially set to about 1 in unit of D/a^ and is decreased by a step of 
AT = 0.0025. At each temperature, states are sampled over r = 5 x 10"^ MC steps. The 
results described below are insensitive to this cooling rate in the sense that values of the 
freezing temperature T*, for example, do not change significantly even when the cooling rate 
is set faster by several times. 

In the present paper we examine a freezing parameter, S, in order to study freezing 
characteristics of the dipolar plate. It is defined as 

s = ^T.m)\u (2) 

i 

where N denotes the number of sites, (...) thermal average over the MC steps r, and [...]r 
average over cooling runs. This parameter represents the frozen degree of spins and is 
important as an order parameter of dipolar systems. In the present work, S is measured 
during cooling processes for the dipolar hexagonal plate whose inter-layer distance / is set 
from /o/3 to Iq, with Iq being the inter-layer distance of the FCC structure. 
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FIG. 2: (a) The temperature dependence of the freezing parameter for each /. (b) The freezing 
parameter as a function of / at T = 0.3. 

II. FREEZING CHARACTERISTICS 

In this section we discuss freezing characteristics based on the freezing parameter whose 
temperature dependence is shown in Fig. |2fa). Data points with each symbol represent the 
values of S for each fixed inter-layer distance, /. At T ~ T*, 5 for each / starts to increases 
rapidly and specific heat (not shown) exhibits a peak. Let us examine the behavior of T* 
in order to clarify trends of the freezing characteristics (see also Fig. Q). At I = Iq, T* is 
nearly equal to 0.65. It steadily decreases as / decreases from Iq. When / reaches nearly at 
O.S/q, T* reaches at its lowest value, which is nearly equal to 0.2. As / further deceases up 
to around 0.7/o, T* slowly increases. After / passes 0.7/o, T* stays nearly constant up to / 
of about 0.5/o where T* starts to increase rapidly. When / is further decreased, T* increases 
continuously. 

From the above observation of Fig. Efa), we can expect at least the following three 
ranges of / in which the freezing characteristics of the system differs significantly; O.S/q ^ ^ 
0.5/o ^ / ^ 0.8/0, and / ^ O.S/q. Actually, at the corresponding temperatures such as at 
T = 0.3, we observe the spin melting and refreezing respectively at / ~ O.8/0 and / ~ O.S/q 
driven by the uniaxial compression as shown in Fig. [2lJ^b). These spin melting and refreezing 
are related to changes of the corresponding ground states as shown in Fig. [T] and discussed 
in detail in the next section. 
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FIG. 3: The ground state for (a) O.S/q ^ ^ (in-plane vortex), (b) 0.51q ^ I ^ O.SZq (convex flow), 
and (c) / ^ 0.5/0 (c) (in-plane-af-FMC). In each schematic picture, spin configurations are shown 
by arrows. The symbols A, B, and C in (c-ii) denote sites on the top, intermediate and bottom 
layers, respectively. 

III. GROUND STATE 

Let us first discuss the spin patterns of the ground state (or the state reached by the 
present analysis at the lowest temperatures) of the plates with / ^ 0.8/o- The one of / = /q is 
shown in Fig. El^a). Spins are confined in the x-y plane and thevform a single vortex. Such 
an in-plane vortex state is often observed in 2d dipolar systems P, 0] . In the present dipolar 
hexagonal plate the vortex state consists of six ferromagnetic domains extended from all 
the edges of the plate. In the cooling process to this ground state the from-Edge-to-Interior 
freezing 0] is observed. Both the energetics of the ground st_ate and the origin of this freezing 



characteristics have been discussed in our previous paper [8|. In a whole range of O.8/0 ~ I 
this vortex state is observed at lowest temperatures. When / becomes close to O.SIq, however, 
the in-plane vortex state are gradually broken from around the center of every hexagonal 
layer, which is indicated by the gray region in Fig. IHl^a). In this region spins possess an 
out-of-plane component. 
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The ground state of the dipolar plate with 0.5/o ^ I ^ O.S/q is found to be a convex- 
flow state schematically shown in Fig. Efb-i). The z-component pattern of spins on the 
middle layer of the plate with / = 0.7/o is shown in Fig. Efb-ii), where a stripe pattern of 
ferromagnetic orders is clearly seen. In fact the absolute magnitudes of these ^-components 
are close to unity. Spins on the first (top) and third (bottom) layers lie almost within the 
layers and align to connect smoothly with the neighboring ferromagnetic orders, yielding a 
convex-flow pattern as a whole. At / ~ O.SZq, the stripe pattern is destabilized to the one 
shown in Fig. El^b-iii). We may call the latter as a randomly-distributed convex-flow state. 
Whether the range of /, where this state appears, is finite or not has not been confirmed, but 
so far it has not been observed at the other end of the convex-flow state, i.e., at / ~ 0.8/o- 

For further compressed systems such as the one with / = /o/3, an in-plane state again 
appears as shown in Fig. E^c). Each on-layer pattern of this state, however, has not the 
vortex structure but the anti-ferromagnetic ferromagnetic-chain (af-FMC) order shown in 
Fig. Efc-i); spins on one lattice axis align ferromagnetically directed to the axis, and this 
ferromagnetic chain order alternates in the direction perpendicular to the axis. It is similar to 
one of the (continuous) ground states of the bulk dipolar system on the 2d square lattice . 
In the present hexagonal plate, the on-layer af-FMC orders align ferromagnetically in the 
2;-direction as shown in Fig. |3fc-ii). Spins on two sides of the plate direct within the side 
plane keeping the af-FMC order with the interior spins, while those on the other four sides 
align to connect continuously between the af-FMC pattern on the top layer and that on 
the bottom layer. The present in-plane-af-FMC ground state is confirmed in the range 
0.33^0 < / < 0.46/o at least. We expect one more change of the ground state at a further 
shorter / than 0.33/o as shown in Fig. ^ because at / = the system is a 2d triangular lattice 
layer whose density of spins is tripled as compared with the original one and so the vortex 
state is naturally expected. 



IV. DISCUSSION 

In our recent work the same dipolar hexagonal systems but with the numbers of 
layers, A'^l, up to 15 and with / = /q, Iq/V^ were studied. We observed the in-plane- vortex 
ground state in systems with / = Iq (/q/v^) for all Al's (A^l < 2), and the convex-flow 
ground state for A^l > 3 in systems with / = Zq/v^- In a bulk system with / = Iq, i.e., the 



FCC lattice, the ferromagnetic order has the global 0(3) symmetry Q|. Therefore, when 
the hexagonal boundary is introduced, there appear the in-plane-vortex ground state which 
gains the boundary- induced anisotropic energy including that from all edges. When the 
FCC structure is compressed in the 2-direction, the magnetization prefers to align to this 
direction, and the competition with the boundary-induced anisotropic energy takes place. 
An important finding of the present work is that the critical distance, Icr, at which the ground 
state changes from the in-plane-vortex to the convex-flow ones, is given as — O.S/q for 
A'^L = 3. More interestingly, the spin melting phenomenon is observed at low temperatures, 
such as T ~ 0.3, as / decreases around 1^, correspondingly to the ground state transition. 
Furthermore, when / is further decreased, the spin refreezing phenomenon is found. It 
corresponds to the ground state change to the in-plane-af-FMC one. The occurrence of the 
latter is certainly related to the fact that the present system with / = /n/2 is just a simple 
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cubic lattice whose one of the bulk ground states is the af-FMC one 
the shape effect are not completely understood yet. 

The results so far described are obtained for the hexagonal dipolar plate by varying I 
but with Nl (= 3) and a size of the plate fixed. Although we have not completed MC 
simulations with systematic change of the other parameters, nor have analyzed effects of the 
boundary roughness, the present results are enough to demonstrate that the dipolar systems 
of a finite size have a vast variety of magnetic phenomena, which are hardly thought of for 
ordinary spin systems only with simple exchange-type interactions. 

The present work is supported by the NAREGI Nanoscience Project of the Ministry of 
Education, Culture, Sports, Science, and Technology. The numerical simulations have been 
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Solid State Physics, the University of Tokyo. 
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